Parton energy loss in a classical strongly coupled QGP 

Kevin Dusling and Ismail Zahed 1 

1 Department of Physics & Astronomy, 
State University of New York, Stony Brook, NY 11794-3800, U.S.A. 

(Dated: December 9, 2009) 

Abstract 

We investigate the energy loss of heavy quarks in the gas, liquid and solid phase of a classical 
quark-gluon plasma (cQGP) using molecular dynamics simulations. The model consists of massive 
quarks and gluons interacting as a classical non-relativistic colored Coulomb gas. We show that the 
electric force decorrelates on a short time scale causing the energy loss to be mostly diffusive and 
Langevindike in the cQGP. We find that the drag coefficient changes with the heavy quark mass, 
while the diffusion constant does not. The fractional collisional energy loss is much larger than 
the leading order estimates from a wQGP (weakly coupled QGP) because of the core repulsion. 
Following recent suggestions, we show how the cQGP results can be translated to the sQGP 
(strongly coupled QGP) results in the T = (1 — 3)T C range. 



1 



I. INTRODUCTION 



One possible signal for the formation of the quark-gluon plasma in heavy ion collisions is 
jet quenching. Partons with a large transverse momentum are created from hard collisions 
between the partons of the nuclei involved in the initial heavy ion collision. These produced 
high pt partons will traverse the collision region (possibly consisting of the quark-gluon 
plasma) within the first few fm/c. Depending upon the properties of the medium traversed 
the parton will lose energy resulting in an experimentally observed jet quenching. 

There have been numerous theoretical calculations of both collisional [l-8] and radiative 



131 ] energy loss of partons in a QCD medium and it is seen that both effects play an 



important role 



14H16I] in the transverse momentum region where RHIC is most sensitive to 



jet quenching. Even with these theoretical results at hand there is still no model which can 
explain the data 21(, i.e. the observed quenching at RHIC energies is stronger then most 
theoretical predictions. 

It should be mentioned that most of the energy loss results were computed in a weak 
coupling expansion which one would expect to converge as long as a s < 1. However, it is 
known from lattice results [3] that at the relevant temperatures probed at RHIC a s reaches 
values of ~ 0.5 and it should be checked how well the perturbative solution converges for 
the values of a s probed in these experiments. 

It is also known that the matter produced at RHIC cannot be weakly coupled but instead 
a good liquid. The evidence for the sQGP [22|, |23J] is large and growing but consists of 
the following two points: 1/ the observed collective flows at RHIC can be explained by 
hydrodynamics showing that the dissipative lengths are very short 2/ Binary bound states 



are seen to exist in lattice simulations above T c and are also predicted 



24| using lattice 



interparticle potentials. We should point out to the reader that there has been arguments 



against bound states as well 



25|. 



Since perturbative methods generally fail in explaining strongly coupled systems other 
approaches have to be adopted. For example, first principle calculations of the sQGP have 
been done using supersymmetric extensions of QCD via the AdS/CFT correspondence. 

The approach taken here, as was first discussed in 2y, |28j, is to model the strongly 
interacting quark and gluon quasiparticles as a classical non-relativistic colored Coulomb 
gas. This model is analyzed using Molecular Dynamics (MD) simulations in which real 



2 



time correlators can be extracted. In 26j decorrelation times, diffusion and viscosity was 
extracted for all phases. It was found that when the results were extrapolated to the sQGP 
a diffusion constant of D ~ 0.1/T and viscosity to entropy density ratio of 77/s ~ 0.3 was 



found. In 



271 ] Charmonium evolution in the sQGP was studied in a related framework. 



In this work we examine the energy loss of heavy quarks propagating through the sQGP. 
In section 2 we shortly review the key ingredients of the cQGP. We show how the structure 
factor can be used to discriminate the gas, liquid and solid phases. We also show that the 
electric force decorrelates on a short time scale in the liquid and solid phase, meaning that 
color probes become rapidly diffusive in the cQGP. In section 3, we numerically assess the 
heavy quark diffusive properties and show that they are amenable to a generic Langevin 
description. In the liquid phase, heavy quarks drag with a drag coefficient that is smaller 
the larger the heavy quark mass. The diffusion constant is independent of the heavy quark 
mass. In section 4, the relative energy loss of heavy quarks for the gas, liquid and solid 
phases are assessed. In section 5, we translate the cQGP results to the sQGP ones in the 
window of temperatures (1 — 3)T C . Our discussions and conclusions are in section 6. Some 
useful units for comparison to the sQGP can be found in Appendix A. A comparison of the 
energy loss with kinetic calculations is in Appendix B. 

II. CLASSICAL CQGP MODEL 

As mentioned in the introduction at temperatures close to Tc quarks and gluons become 
quasiparticles with masses on the order of 3T. We can model the sQGP as a system of mas- 
sive non-relativistic particles interacting through longitudinal color electric fields. Magnetic 
effects are suppressed in the non-relativistic limit. The specific Hamiltonian used in our 
model of the sQGP is 



H 



^ 2m, ^ 



QaiQ^j 



+ v a 



(1) 



where 



V„ 



(2) 
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with i,j = l..N a being a sum over all particles of specie a,j3= q,q,g having respective 
particle number N a . 

The first term in the above Hamiltonian is the standard kinetic energy term. The second 
term is the colored coulomb interaction. Again we have neglected any chromomagnetic 
interaction in the non-relativistic limit. Non-perturbative effects due to magnetic charges 
are discussed in 



29 



34j. A short range repulsive potential was added by hand and is needed 



to give stability to the simulation and can be argued to mimic the effect of a quantum 



localization energy. A more detailed study o 



classical MD simulations can be found in 



35]. 



the quantum corrections to potentials used in 



The equations of motion can be derived from the usual Poisson brackets (O a i = {H, O a i}) 
where O ai is the phase space coordinate of either position (x ai ), momentum (p ai ) or color 
(Qai) of particle ojj. 

The strength of the interparticle interaction is classified in the context of traditional 
electromagnetic plasmas, using the dimensionless parameter T, the ratio of the potential to 
kinetic energy: 



(3) 

a ws T 

where Ze, aws, T are respectively the ion charge, the Wigner-Seitz radius aws = (3/47rn) 1 / 3 
and the temperature. One usually defines the weakly coupled or gas regime for T < 1, a 
liquid regime for r ~ 1 — 10 and a strongly coupled or solid regime for for T > 10. 

The gas, liquid and solid nature of the cQGP can be seen by looking at the structure 
factor for the different phases. In Fig. [1] we show the distinct-correlation function (Gd) 
defined as: 



G d (x,t) = j^< J2S{x + Xi{Q)-£j{t)) > (4) 

for a gas, liquid and solid phase having T rs 0.1,3,160 respectively. Gd characterizes the 
probability to find two distinct particles having a separation r at a time t. One can see the 
that correlation among particles increases as one goes from the gas to liquid to glass phases. 
It turns out, as discussed in appendix A, that we expect the properties of the sQGP to be 
consistent with the cQGP having r ~ 3. 
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FIG. 1. (Color online) correlation function for V ~ 0.1,3.0, 160 and t = 0. 

The nature of the decorrelation times can be numerically investigated. The color electric 
forces decorrelate on a short time scale in comparison to the velocity decorrelation, meaning 
that color probes whether heavy or light are readily diffusive in the cQGP whether gas, liquid 
or solid. The electric force follows from (CQ) as Qi a F ia = —dH/dxi a , and the decorrelation 
function 
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G ia {t) =< F ia {t)F ia {0) > I < F ia (0) 2 > (5) 

is shown in Fig. [2] for the gas, liquid and solid phase. The color electric force decorrelates 
rapidly in the gas and liquid phase, and more slowly in the ordered solid phase. In the liquid 
phase with T = 3 the decorrelation time is tp = r/2 in simulation units. We translate these 
units to physical units in Appendix A. In particular tp ~ 1/10T for the sQGP. 



III. DIFFUSING JETS 



To analyze the evolution of jets in the cQGP at strong coupling, we will use molecular 
dynamics simulations (see {26] for more details) of the cQGP by integrating the above 
equations of motion of 64 particles confined to a box with periodic boundary conditions 
and mirror cubes on all sides. When equilibration of the system is reached an external 
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FIG. 2. Color electric force-force decorrelator for T ~ 0.1, 3.0, 160 



particle is added to the simulation with a given momentum and mass (M) greater than 
the quark and gluon quasiparticle mass (m). The evolution of this probe particle is tracked 
throughout the evolution of the system and it phase space coordinates are recorded. This 
procedure is repeated until enough measurements have been made to make a statistically 
significant distribution of any phase space quantity. For a given momentum bin we use 
m 1500 measurements. Statistical errors are given by yN of the number of runs. 

In Fig. [3] we show two examples of the evolution of the probe particle's momentum 
distribution function, f(p), as measured in the MD simulation in the liquid phase with T = 3. 
The first shows the evolution for a heavy quark mass with M = 2m and pi n a = 10 (large 
initial peak) and the second for a heavy quark mass with M = 10m and Pi n u = 20 (large 
initial peak). Each figure shows the initial momentum distribution and the distribution 
at two later times. As the time of flight or distance increases a decrease in the average 
momentum of the probe as well as a broadening of the distribution function is seen. 

Since the color electric forces decorrelate promptly as shown in Fig. [2J the massive jet 
enters rapidly a diffusive regime. Assuming homogeneity in space, the diffusion is mostly in 
momentum space. This is amenable to a non-relativistic Fokker-Planck equation (Ornstein- 
Uhlenbeck type) whereby the lowest two moments, i.e. drag and diffusion, are dominant. 
Restricting our analysis to the one-dimensional case for simplicity [stJ 
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FIG. 3. (Color online) Momentum distribution function, f(p), shown for various values of time 
spent in the liquid cQGP. Left: Pmit = 10 and M = 2m, Right: Pmit = 20 and M = 10m. 
All quantities are expressed in simulation units and measured from a simulation with T = 3 
corresponding to a temperature T « 7. 



% = ^m + ~wpm (6) 

where the lowest two moments M\ and Mi are related to the drag and diffusion constant 
respectively, 

Mi(p) = ^ = ~V(P)P 

M 2 (p) = { -^p- = K L (p) (7) 

These two moments can be measured directly in our MD simulation. The first moment, 
Mi, is related to the average momentum loss where r)(p) is the drag coefficient and as is seen 
in fig. H] is independent of p for initial momentum greater than the thermal momentum of 
the system. We therefore take 77 (p) = 0.35,0.2 and 0.1 (in simulation units) for the various 
cases of the heavy quark mass having M = 2m, 6m, and 10m respectively. Noticeably, the 
drag coefficient 77 (p) decreases as the jet mass increases. The second moment, M2, is related 
to the longitudinal momentum fluctuations (kl) and is oc E as shown in fig. HJ Within 
uncertainties k l is independent of the heavy probe mass and we take k l = 0.25E. 



Using the parameterizations of r\ and kl defined above the evolution equation [6] can be 
solved numerically as shown in fig. [3] as solid curves. Within the uncertainties of both the un- 
derlying model (soft collisions, one-dimensional) as well as our parameterization of transport 
coefficients the Fokker-Planck analysis shows good agreement with the MD distributions. 
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FIG. 4. (Color online) Measured transport coefficients r/ and kl for three cases of the probe 
particle's mass. All quantities are expressed in simulation units and measured from a simulation 
with r = 3. 



IV. ENERGY LOSS 

From the momentum distribution function we can evaluate the mean energy of the ex- 
ternal parton after traversing a distance L as: 



(E}= / E f(p, L) dp (8) 
Jo 

and define the fractional energy loss as: 



AE = E - (E) 
E Eq 

In figure [5] we show the fractional energy loss as a function of length (in simulation units) 
for three different coupling parameters. In all cases the external particle has a mass ten 

8 




2 4 6 8 10 5 10 15 20 25 30 35 40 



L E 

FIG. 5. Fractional energy loss as a function of length and energy measured from simulations having 
coupling parameters of T = 0.8, 3.0, 30. All quantities are expressed in simulation units. 

times that of the quasiparticle mass (i.e. M = 10 in simulation units). For all three cases 
an approximately linear rise in energy loss as a function of length is seen which is expected 
for collisional loss only. We also show the fractional energy loss as a function of initial 
energy (in simulation units) for three different values of coupling parameters. As expected, 
at thermal energies, the fractional energy loss is negative because a particle with v — 
can only gain momentum in collisions. One sees that the fractional energy loss remains 
constant at high enough energies. The fact that the energy loss is proportional to E at high 
enough energy instead of increasing logarithmically as in the case for Coulomb collisions is 
due to the core potential (eq. [2]). Since the different coupling parameters T are modified by 
changing the system temperature we don't expect to see differences in the energy loss when 
the probe's momentum is much greater than the thermal momentum p ^> \JlMT of the 
quasiparticles. 



V. COMPARISON TO SQGP 



In order to compare the results from 
at RHIC we follow the arguments of 26 



;he MD simulation in the cQGP to the sQGP 



281 ] as summarized in Appendix A. Specifically, 



we consider a plasma at a temperature T = 250 MeV~ 1.5T^ corresponding to a coupling 
parameter T = 3. The unit of mass is rescaled using [m] = [3T] « 0.75 GeV, so a heavy probe 
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with mass 2, 6 and 10 times the quasiparticle mass as used in the simulation corresponds 
to masses of 1.5 GeV, 4.5 GeV and 7.5 GeV respectively in the sQGP. As mentioned earlier 
the plasma consists of quasiparticles with m ~ 3T pa 0.75 GeV at T = l.hTc- The drag 
coefficient has units [77] = [1/r] = [5. IT] pa 6.5 fm -1 . Assuming that the drag is independent 
of p as we showed for a large range of momenta we find that rj = 2.3, 1.3, 0.65 fm -1 for probe 
masses of M = 1.5,4.5, 7.5 GeV respectively. 

It is also useful to look at the energy loss after the probe particle travels a finite length L 
(here taken to be about the length of our simulation region or 4.5 simulation units) which 
corresponds to 1.2 fm in the sQGP. In fig. [6] we show the fractional energy loss as a function 
of the probe particle's initial momentum as calculated from the MD simulation (points with 
error bars). For comparison we also show the leading order in a s collisional energy loss of 



171 ] . The sensitivity to 



the scattering of heavy quarks off of massless quarks and gluons 
the schematic form of the core potential is smaller at larger values of momentum. This is 
verified by a comparison of the MD simulation results to analytic calculations preformed in 



^ = _s™^ (1 + rv )(1 _i^ ln i±^ )ln£ ^ (10) 

dx 3v v 6 ;v 2v l-v 1 q mn v ; 



where the lower cutoff of the momentum transfer is taken as the Debye mass: q mm ~ 2T and 
the upper cutoff q rnax ~ y/ ATE is taken from [2]. In order to compare the above equation 
to our simulation results we take AE/E pa L/E ■ (—dE/dx) where L = 1.2 fm. 



VI. DISCUSSION AND CONCLUSIONS 

First, it should be mentioned, that the energy loss in equation [10] acts as an upper limit 
for other energy loss calculations as presented in Comparing the MD simulation 

results to the leading order results in eq. [10] one immediately sees a much larger fractional 
energy loss for the sQGP. As much as a factor of five in the high momentum region. The 
main physical differences that account for the large increase are: 1/ the use of quasiparticle 
quark and gluon masses (m ~ 3T) much greater than the current quark masses 2/ the use 
of the hard core potential to mock up the effects of quantum repulsion at short distances. 
3/ any non-local many body interactions induced by the coulomb and core potential are 
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FIG. 6. Fractional collisional energy loss as a function of momentum for charm, bottom and heavy 
quarks (top, middle and bottom figures respectively) in the sQGP. 
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resummed to all orders using MD. 

One should note that if the current lattice results atTw l-5Tc are modified the same 
MD simulation results can be used by simply applying a different rescaling of the units 
consistent with the data. Since the fractional energy loss is mostly constant as a function 
of momentum, a rescaling of the momentum will not change the results in the high energy 
region. However, both the length of the medium and the mass of the heavy probe particle 
will need to be adjusted accordingly. The change in length unit will also modify the result 
for the drag coefficient. 

In conclusion we have computed the parton energy loss in a strongly coupled quark-gluon 
plasma using a classical molecular dynamic simulation. The model consists of massive 
strongly interacting quark and gluon quasiparticles at T = (1 — 3)T C interacting as a classic 
non-relativistic colored Coulomb gas. We find that the fractional collisional energy loss is 
much larger compared to the leading order estimates in a wQGP. 
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Appendix A: Units and Comparison to sQGP 



The equations of motion provided above are integrated over time given some initial con- 
figuration of phase space coordinates. It is convenient to run the evolution in simulation 

]. First, the unit of length is set by the minimum 



units which we now discuss following 
of the potential which has the form: 



26 



V 



9_ 

A 



Q-Q 



A 1 ,A, 



(Al) 



where A = r min (for any value of d) sets the basic length scale in which all distances are 
measured. The time unit is set by the plasma frequency of the system: 
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and the unit of mass is denned as the particle mass. For example, in a simulation of only 
one particle specie all masses are equal to one. 

In order to extract results from the cQGP simulation about the physical sQGP, all that 
is required is a re-scaling of the three basic units of length, time and mass. For now this 
is done at a temperature of T = 1.5 — 3T C where most lattice data is available. The unit 
of mass is set by the mass of the quark and gluon quasiparticles taken from lattice data: 
m ~ 3T at T = 1.5T C as discussed earlier. 

The effective interparticle potential is given by: 



2mr 2 r (^) 

where C is the pertinent Casimir for quarks and gluons. The length unit in the cQGP is 
set by the minimum of the potential, A = h 2 /mCa s . With a s ~ 0.5 at these distances 
and averaging over the color casimir for quarks and gluons assuming that all three species 
(g, q, q) are equally represented, we have < a s C >~ 1 within uncertainties of the model. 
This leads to A ~ ^ in units where h = c = 1. 

Finally the time unit is given by { 4nn< ^ sC> ) ■ The density of quasiparticles (n) is 
estimated as the density of black body radiation photons multiplied by the effective degrees 
of freedom: n ~ (0.244T 3 )(8 + 6Nf) rs 6.3T 3 . Then the time unit is measured in r ~ ^4^. 



Appendix B: Comparison with Kinetics 

In order to understand the collisional energy loss better we compare the results from 
the full MD simulation with a simple kinetic calculation. We assume that the heavy quark 
undergoes two body scattering with a plasma of massive quasiparticles of density n. When 



the heavy particle passes through the p 



parameters. The energy loss is given by 39(: 



asma it sees quasiparticles at all possible impact 



^— = 2im / T{b)bdb (Bl) 
ax Jq 

where T(b) is the energy transfer to the heavy particle in a two body collision with a 
quasiparticle at rest at an impact parameter b. Fig. [7] shows the energy transfer as a function 
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FIG. 7. Energy transfer as a function of impact parameter for a heavy quark with Pj n = 15 and 
mass M=2m, 6m, and 10m 

of impact parameter for three heavy quark masses and an initial momentum of P« n = 15 in 
simulation units. The results were calculated for the same coulomb+hard core potential as 
used in the MD simulations. For comparison the solid curve shows the analytic result for 
scattering from a coulomb potential in the limit that m/M — > 0: 

2cy 2 1 

*W - -TjgfTv (B2) 

For distances larger than one T(b) increases towards the Coulomb the mass is in- 

creased as is expected. For distances smaller then one there is a large enhancement in the 
energy transfer due to the strong repulsive core of the MD potential. 

In fig. [8] we show the drag given by equation IB II compared to the MD simulation results 
already presented from fig. HI The comparison is only shown for velocities larger then the 
thermal velocity where equation IB1I holds. There is good agreement between this simple 
model and the full MD results. 
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